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Abstract: We construct a simple holographic QCD model to study nuclear matter to 
strange matter transition. The interaction of dense medium and hadrons is taken care of by 
imposing the force balancing condition for stable D4/D6/D6 configuration. By considering 
the intermediate and light flavor branes interacting with baryon vertex homogeneously 
distributed along R s space and requesting the energy minimization, we find that there is 
a well defined transition density as a function of current quark mass. We also find that as 
density goes up very high, intermediate (or heavy) and light quarks populate equally as 
expected from the Pauli principle. In this sense, the effect of the Pauli principle is realized 
as dynamics of D-branes. 
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1. Introduction 

One of the challenging current problems in hadron physics is to elucidate the behavior 
of dense matter under extremely high-density environments, for a review, e.g. see It 
is theoretically expected that, at very high baryon densities (even at low temperatures), 
chiral symmetry is likely to be restored, and that baryon matter can be converted into 
quark matter. Various studies also suggest the possible formation of a kaon condensate 
at high densities. The existence of quark matter and/or a kaon condensate can have 
important consequences for the structure of compact stars and for the cooling behavior 
of a remnant star after supernova explosion and the subsequent formation of a neutron 
star. Thorough understanding on dense nuclear matter is also important in the physics of 
relativistic heavy- ion collisions. 

In this study, we focus on a specific aspect of dense matter: transition from nuclear 
matter to strange matter. This transition is essential to understand relevance of kaon 
condensation in neutron star. This is because the presence of strangeness matter tends to 
hinder the formation of kaon condensation basically due to the Pauli exclusion principle. 
Here kaon condensation means K~ condensation. Since K~ is composed of u and s, strange 
matter will expel the strange quark and so the kaon condensation. Moreover, according 
to previous studies, the critical baryon number density for the nuclear matter to strange 
matter and that for the onset of kaon condensation are not very different from each other, 
~ (2 — 4)po; where po is the normal nuclear matter density ~ 0.17 fm~ 3 ~ In 
conventional approaches, however, when one estimates the transition to strange matter, 
a bit large uncertainty comes in due to lack of robust information on hyperon coupling 
constants, see [|[ for a review. In this study, we take a first step towards this direction in 
holographic QCD. Recent developments based on AdS/CFT ||, |||, ||] renders a new tool 
to study dense matter in the framework of a holographic model of QCD || 0], see || for 
a review. There have been many studies on dense nuclear matter |l(J in the holographic 
QCD. 
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To study the transition we introduce two flavor D6 branes which correspond to light 
(u or d) and intermediate (strange) quarks respectively and spherical D4 brane with Nc 
fundamental strings. The fundamental strings can be attached on a light quark D6 brane 
and/or intermediate mass quark (strange) D6 brane. By solving DBI action and applying 
the force balance condition, we find stable configuration of D4/D6/D6 system. After 
finding minimum energy configuration, we calculate the ratio of up and strange quarks in 
the system function of baryon density. 



2. Baryon vertex 

In this section we discuss baryon vertex (spherical D4 brane with Nc fundamental strings) 
in confining background following [|ll]] 

The non-supersymmetric geometry for confining background of D4 in Euclidean sig- 
nature is given by 

ds 2 = (|) ^ (^dxv + f(u)dxl) + (^) 3/2 (J^ + u 2 dn 2 



[A 3/4 2irN c /U K 



6<P = 9s { 1 ) ' Fi = ~7^ e4 ' = 1 " ' R3 = ^sNJl (2.1) 

This background is related to the black hole solution of D4 brane by the double Wick 
rotation. The Kaluza-Klein mass is defined as inverse radius of the X4 direction: M^k = 

I f ^ . While Ukk,9s, and R are bulk parameters, Mmm and g\ M are the parameters of 
the gauge theory. These are related by 

U K K = l\M KK ll # 3 = -^-, \ = g ^ M N c . (2.2) 



ys 2-kI s N c Mkk 9 ** s ' 2M KK 

Introducing a dimensionless coordinate £; ^p- = j0jwp, , we obtain the background geome- 
try 

ds2 =(X) (dt 2 + dx* + f(U)dxl) + (^J (£j (df + fdnl). (2.3) 

Here U and £ are related by {m^\^ = \ {^ /2 + ^72 ) • 

A baryon in three-dimensional theory corresponds to the D4 brane wrapping S 4 on 
which N c fundamental strings terminate. In this configuration, the background four-form 
field strength couples to the world volume gauge field via Chern-Simons term. 

The background metric ( |2.3| ) can be written as 

ds 2 1 



R 



\ (dt 2 + fdx 2 + dx 2 ) + R 3 / 2 Vu(^ + d8 2 + sin 2 0dfi§V (2.4) 



We take (t, 6 a ) as a world volume coordinate of D4 brane and turn on the U(l) gauge field 
on it, F t Q ^ 0. For simplicity, we assume that the position of D4 brane and the gauge 
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field depends only on 9 i.e. £ = A = At(0), where 9 is the polar angle in spherical 
coordinates. The induced metric on the compact D4 brane is 



ds 2 



D4 



TI\ 3 / 2 

dt 2 +RV 2 VU 



1 + T2 ) d ° 2 + Sin2 ^3 



(2.5) 



where £' = d^/d9. The DBI action for single D4 brane with iV c fundamental strings is 
given by |§ 

Sda = — A*4 y e~*\/det(5 + 2vrQ'F) + /x 4 ^ A G (4) 
= r 4 / dtdd sin 3 -v/w+ 3 (£ 2 + £ /2 ) - F 2 + 3i< 



dtCo4, 



(2.6) 



where 



r 4 = /M^&j # 
~ 2vra'Fi 2 4 / 6 



i -t&Ukk _ N c Ukk 



Ukk 



2 2/3 2 8 / 3 (2vrZ 2 ) 
2 2 / 3 

A = • 27ra A 



(2.7) 



with (j± = 1 ± £ 3 . The dimensionless displacement is defined as follows; 



dC 



Dl 



sin 3 OF 



dF 



V/<4 /3 )(£ 2 +£' 2 )-^ 2 

= -D(9). 

Then the equation of motion for the gauge field is 

d B D(6) = -3sin 3 #. 
By integrating above equation, we get 

1 Q 

D(6) = 2(2i/ - 1) + 3(cos B - - cos 3 9). 

3 



(2.9) 



(2.10) 



where the integration constant v determines the number of fundamental sting; vN c strings 
are attached at south pole and (1 — v)N c strings at north pole. 
After Legendre transformation, we obtain 'Hamiltonian' as 



~ dC Di 

riDi — r — — ~ '-'DA 

oF 



r 4 J d9^u; A f(¥ + i' 2 )^D(9) 2 + sin 6 9, 



(2.11) 



We solve the equation of motion for 7i £> 4 numerically. We set v = since we assume that all 
fundamental strings are attached at north pole 9 = n. Then we impose smooth boundary 
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Figure 1: Shape of D4 brane for different Re- 



condition £'(0) = and £(0) = at 8 = 0. Numerical solutions are parameterized by 
initial value of £o- The solutions corresponding to different £;o's are drawn in Fig. |l|. 

If we denote the position of the cusp of D4 brane by U c , the force at the cusp due to 
the D4 brane tension can be obtained by varying the Hamiltonian of D4 brane with respect 
to U c while keeping other variables fixed; 



&H 



dU c 



N C T F s % - sc (2.12) 



fix other values 



where Tp = £ is tension of fundamental string. The tension at the cusp of D4 brane 
is always smaller than the tension of the Nq fundamental strings. Therefore, if there are 
no other object, the cusp should be pulled up to infinity and the final configuration of D4 
brane would be 'tube-like' shape as in p 



3. Holographic transition to strange matter 

We begin with an simplified description of the nuclear matter to strange matter transition. 
Figure ^ shows a simple view of the transition. Here we set aside issues like charge neutrality 
and /3-equilibrium of the dense matter and consider two light quarks (u, d quarks) and one 
intermediate mass quark (s quark). The vertical axis of Fig.|| is roughly the number of 
quarks in the ground state since n q ~ kp, where n q is the quark number density. In low- 
density regime, we would have only u and d in our system since the mass of the strange 
quark, m s , is roughly a few ten times bigger than that of light quarks, m s /m u ~ 50, 
m s /md ~ 20, and so cost too much energy to be piled up in the ground state. As we increase 
the number density, the chemical potential of light quarks (Hu,d) become comparable with 
the mass of the strange quark, and then system could lower its ground state energy by piling 
up some strange quarks if [i u ^ > m s . This argument goes also with baryons. Instead of 
transition from u, d quark matter to u, d, and s matter, it could be a transition from 
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nuclear matter composed of nucleoli to strange matter of nucleon and hyperons like A. As 
it stands, this is too simple. To be realistic we have to include interaction energy, charge 
neutrality of the matter, /^-equilibrium, etc. 



u d 



d 



(a) (b) 
Figure 2: Schematic picture for the nuclear matter (a) to strange matter (b) transition. 



Now, we delve into the transition in terms of a holographic QCD. For simplicity, we 
will again ignore the charge neutrality and /3-equilibrium of the matter at hand, relegating 
them to a future study. Now we consider the system with two flavors, one light and one 
intermediate mass quarks. To introduce two flavors, we put two probe D6 branes in the 
background. The the bulk metric (|2.3|) can be written as 



3 /2 3/2 2 

ds 2 = (|j (dt* + dx* + f(U)dxt) + ($\ (J) (dp' + p'dnl + dy* + y 2 # 2 ) , 

(3.1) 

where D6 brane world volume coordinates are (t,x, p,8 a ). The embedding ansatz is that 
only y depends on p (we set = 0). The induced metric on a single D6 brane is 



ds 



D6 



|) ((ft 2 + dx') + (|) 2 [(1 + y W + p 2 dnl] , (3.2) 



where y = dy/dp. 

We also introduce U (1) gauge field on D6 brane, which is coupled to the string end 
point. The DBI action for the single D6 brane is 



Sdg = J dtC D6 = -p 6 J e-VdetG? + 27ra'F) 



where 



-r 6 J dtd PP 2 ujf^f(l + y2)-F2, 



r 6 = >0 25 -^, F^-^f-. 



(3.3) 



(3.4) 



We define dimensionless quantity Q from the equation of motion for F; 

dS D6 _ p 2 Jj 3 F 



dF 



4 /3 (l + y 2 )-F 2 



Q. 



(3.5) 
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The Hamiltonian is connected to the number of point sources (number of fundamental 
strings) Q by 

UkkQ 



Q 2-2 2 / 3 7ra'T 6 ' 
The Hamiltonian can be obtained by the Legendre transformation; 



(3.6) 



T~(-D6 



OF 



T6 



r 6 J d P v(p)Vi + y 2 



(3.7) 



To solve the equation of motion, we impose appropriate initial condition to D6 brane. 
We are considering two D6 branes that are connected to a D4 brane with fundamental 
strings. As discussed in [11], the tension of fundamental strings is always larger than that 
of D-branes. Therefore, two D6 branes are pulled down and spherical D4 brane pulled up 
until length of fundamental strings becomes zero. Finally, the position of the cusp of D6 
branes should be located at the same position of the cusp of D4 brane £ c . We consider 
Qi fundamental strings attached on one D6 brane and Q 2 strings attached on another D6 
brane. We also denote the slope at cusp of each brane as g/c and yi 2 ^ . The force at the 
cusp of D6 branes can be obtained as 



dH(Qi) m 



dU, 



+ 



dH(Q 2 ) D6 



•X 1 ) 



2ira' V 1 - 



+ 



Q2 ri+i 



i + ^ 1)2 



-3 



2ttq' V 1 - £, 



—3 



•,(2) 



1 , -(2)2 
1 + Vc 



(3.8) 



To make the whole system stable, the force at the cusp of D4 brane should be balanced to 
force of D6 branes; 



^f D4 = f«(q 1 ) + f( ) 2 6 ) (q 2 ) 



(3.9) 



where Q = Q\ + Q 



2- 



.-.(0 



Rewriting Qi and ijc by using new parameters a and f3, 



Qi = (l-a)Q 



y, 



(i) 



Q 2 = aQ 2 , 



(3.10) 



the force balance condition (|3.9|) becomes 



& (1 - 



+ 



a(3y, 



(1)2 



1 i «2-(!)2 

l + p l yl 



(3.11) 



We show solutions that satisfy the force balance condition above in Figure.|3](a). 
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2=8 




Here we fix asymptotic values of D6 branes which are correspond to quark mass to 
be mi = 0.1 and m<i = 5. We call these branes as up and strange quarks brane for 
convenience. We note here that mi is dimensionless and it is related to the real quark mass 
by the following relation 

rn.XMKK (312) 
y7r 

We will use this relation to discuss the parameter fixing. 

The behavior of solutions change depending on the density and £q. If we fix £o> then 
from the equation of motion of D4 brane, the cusp point ^ c and slope at cusp £' c are 



automatically determined. However, from eq. (4.1), for given ^ c and slopes of two D6 
branes which give mi = 0.1, ni2 = 5 cannot be uniquely determined. In fact, there are 
infinite set of a and (3. So we have to choose one embedding by minimum energy condition. 
The total energy of this system can be written as follows; 



Etot — -Q-HDi + Ti-DeiQi) + 7~Ld&{Q2) 



T6 



+ £ 6 (Ql) + E & {Q 2 ) 



(3.13) 



where Ei is numerical integration of each 'Hamiltonian' in ( p. 11 ) and ( |3.7D without overall 



constant T4 and tq. The a dependence of total energy is drawn in Fig. |3](b). In this figure, 
we can see that the total energy of system is minimized at at a = 0.361. In other words, for 
Q = 8 case, 36.1% of total quarks of ground state (or Fermi sea) are occupied by strange 
(intermediate mass) quarks. For each Q, we find corresponding a by imposing the energy 
minimum condition. 

For small value of Q, a dependence of total energy have different behavior Fig. ||](b) as 
shown in Figure. f|(a). In this figure, total energy monotonically decreases as a decreases. 
Finally, the minimum energy configuration is at a = 0, it means that at small density, 
strange quark cannot come into the system. The embedding in this case is drawn in 
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Figure. ||. As we increase the density, at a certain density, a = is not a minimum energy 
embedding anymore. From Q ~ 4.2, non-zero a embedding has minimum energy, and the 
value increases as density increases. We also check the large density behavior of embedding: 
when Q is large, the ratio of strange quark to up quark seems to go to 0.5, for example 
when Q = 500 the value of a is around 0.49. The final result is drawn in Figure. |f|(a), 
which shows sharp transition from nuclear matter (a = 0) to strange matter (a ^ 0). It 
is interesting to note that similar tendency has observed in QCD-rooted model studies, for 



instance see [13]. 





(a) 



(b) 



Figure 5: (a) Density dependence of a, the fraction of the strange quarks, (b) Density dependence 
of ro, a measure of the baryon mass, for mi = 0.1 and to 2 = 5. 



Next, we calculate the value of r$ for each embedding which is proportional to the 
energy of spherical D4 brane. As we discussed in [11|, this value can be interpreted as 
mass of baryon. For small value of Q, the behavior of r$ is the same as [|ll]] because in this 
region only one probe brane touch the baryon vertex - as Q increase, po decreases first and 
then increases. After Q ~ 4.2, two probe branes are attached to baryon vertex and ro is 
increase. The behavior of ro as a function of Q is drawn in Figure. ||(b). 
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Intuitively as mj/mi increases, the transition density Q c from zero a to finite a, or 
from nuclear matter to strange matter, should increase. We check if our D4/D6/D6 system 
follows this expectation. For this we used three different values for 771-2, Tn,2 = 2,3,5, with 
mi = 0.1. The results are summarized in Fig. || to show that our D4/D6/D6 system 
respects the intuition. 



mi =0.1 




Figure 6: (a) Density dependence of a with different m^/mi. All lines saturate a = 0.5 for large 
Q. (b) p c vs. m 2 /mi. 



Though our D4/D6/D6 system is surely far from a realistic dense QCD matter at this 
stage, we convert the transition Q c obtained in this study into a baryon number density in 



QCD. From (}L2[ ) and Q3.4J) , we can calculate the baryon number density in terms of Q as 
follows; 

= Nb _ _Q_ 
9 ~ V 3 N C V 3 
2 • 2 2 / 3 



For mi = 0.1 and m-2 = 5, we obtained Q c ~ 4.2. Empirically the transition baryon 
number density is ~ (2 — 4)po> where po is the normal nuclear density, from various model 
studies based on QCD. To compare our result with the empirical transition density, we 
choose A = 1.2 and Mkk = 1-5 GeV. With this choice we obtain p c ~ 1.96/3o for Q c ~ 4.2. 
From ( p. 12 ), corresponding quark masses are M\ ~ 6.4 MeV and M2 ~ 318 MeV, which 



are close to the mass of up and strange quarks in QCD. At the end of the day, however, 
on completion of our D4/D6/D6 model for realistic dense matter, we may have to fix the 
values of A and Mkk ah initio by considering fluctuations of D6 branes and by doing meson 
spectroscopy. 

So far we consider mi = 0.1 and m<i = 2,3,5. For fixed ratio between mi and m-2, 
however, we can choose several different embeddings. 

• Case 1: mi << Ukk = 1 < m2- This is the case we discussed above. 

• Case 2: both of m q are smaller than Ukk, such as mi = 0.01 and m-2 = 0.5. In 
this case, brane embedding is drawn in Figure. 0(a). The density dependence of a is 
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drawn in Figure. 0(b). In this case, the value of a seems to saturate to 0.5 for any 
non-zero density. That is, it does not go to zero even in extremely small density. We 
can understand this behavior in geometrically. If both of m q is smaller than Ukk, 
the difference of geometry detected by each brane is very small. So, the each brane 
shares nearly equal number of quarks. 

• Case 3: both of m q are larger than Ukk- Here we expect similar behavior as in case 
2. Two quarks will populate evenly. 

Finally we remark that in our model the value of a depends on not only the ratio of mi 
and mi-, but also each value of mi and m,2- 



0.50 - 
0.45 - 
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(a) (b) 
Figure 7: (a) Brane embedding for m\ = 0.01 and m-i — 0.5. (b)Density dependence of a. 



4. Summary and Discussion 

In this paper, we constructed a simple holographic QCD model to study nuclear matter 
to strange matter transition. The interaction of dense medium and hadrons was taken 
care of by imposing the force balancing condition for stable D4/D6/D6 configuration. We 
considered the oversimplified model where only one intermediate and one light flavor branes 
are interacting with baryon vertex which is homogeneously distributed along i? 3 space. We 
imposed the energy minimization condition and found that there is a well defined transition 
density as a function of current quark mass and that the transition density increases as the 
ratio mi/fni increases. We also showed that at very high densities, intermediate (or heavy) 
and light quarks populate equally as expected from the Pauli principle, (add a few lines 
and lower the voice) So, we may conclude that in our study the effect of Pauli principle is 
realized as dynamics of D-branes. 

We could lessen the oversimplification a little bit by considering three flavors with 
m u = rrid < m s , which could be realized by considering the light quark favor brane to have 
weight 2 relative to the intermediate one in the force balancing condition. That is, if we 
change it to 

^-F DA = 2F§l(Q 1 ) + F%l(Q 2 ) . (4.1) 
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In this case, it is expected that the asymptotic value of a would be 1/3 instead of 1/2. To 
confirm this we take Q = 500 and obtain a ~ 1/3. 

It will be interesting to consider fluctuations in our D4/D6/D6 system to study meson 
masses in iso-spin asymmetric matter, which will be reported elsewhere. 
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